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STRONGLY  TURBULENT  STABILIZATION  OF  ELECTRON 
BEAM-PLASMA  INTERACTIONS 


I.  Introduction 

The  interaction  of  relativistic  and  non-relativistic  electron 
beams  with  a  plasma  remains  an  important  topic  of  research.  In 
addition  to  the  usual  applications  both  to  space  plasmas  and  the 
heating  of  linear  devices  in  the  laboratory,  a  resurgence  of  interest 
in  this  problem  has  occurred  due  to  possible  applications  in  toroidal 
magnetic  confinement  devices.  For  example,  Mohri  et  al . ^  have  injected 
relativistic  electron  beams  into  toroidal  devices,  and  Papadopoulos  et  al.2 
have  proposed  the  in  situ  formation  of  relativistic  electron  beams  in 
tokamaks  as  a  current  driver,  for  long  pulse  operation,  and  for  supplemen¬ 
tary  heating.  Finally,  it  should  be  mentioned  that  several  of  the  spheromak 
concepts  require  relativistic  electron  beam  drivers.  Any  assessment  of 
these  concepts  requires  a  good  theoretical  basis  of  prediction  of  the 
beam  relaxation  process.  It  is  our  purpose  in  this  work  to  provide  such 
an  understanding  based  upon  computational  means. 

The  failures  of  the  quasilinear  (or  any  weak  turbulence  theory)  to 

describe  the  beam  relaxation  process  is  well-known  in  the  literature3-3, 

and  the  strong  turbulence  theory  has  been  successful  in  resolving  many 

of  the  discrepancies  between  theory  and  experiment.3-3  However,  in 

view  of  the  complexity  of  the  problem,  several  ad  hoc  simplifications 

are  commonly  made  in  the  name  of  analytic  tractability .  We  present, 

here,  the  results  of  a  numerical  solution  of  the  strong  turbulence 

equations  which  describe  the  temporal  evolution  of  a  beam  streaming 

through  a  plasma.  The  model  w;  employ  is  one-dimensional,  which  is 

easily  justified  in  the  presence  of  even  a  modest  magnetic  field.  We 

neglect  the  self-consistent  interaction  of  the  waves  with  the  beam 
Note:  Manuscript  submitted  October  9,  1979. 
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particles  by  the  assumption  of  a  constant  linear  growth  rate  for  the 
beam-plasma  instability.  This  assumption  can  be  justified  a  posteriori , 
and  will  be  discussed  in  more  detail  in  the  concluding  section. 

Following  the  presentation  of  the  simulation  results,  we  develop  a 
phenomenological  model  to  describe  the  beam  relaxation-process  in  an 
attempt  to  reproduce  the  scalings  of  the  fluctuation  quantities  with 
the  linear  growth  rate. 

9 

In  addition,  Manheimer  and  Papadopoulos  have  demonstrated  that 
the  equations  which  govern  the  nonlinear  coupling  between  Langmuir  and 
ion  acoustic  waves  are  formally  equivalent  to  Zakharov's  equations^  for 
Langmuir  solitons.  As  a  consequence,  it  is  reasonable  to  expect  the 
formation  of  soliton-like  structures  to  occur  by  means  of  this  mechanism, 
and  such  a  ’esult  is,  in  fact,  clearly  shown  by  the  numerical  solution 
of  the  dynamical  equations. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we 
discuss  the  numerical  simulation.  The  basic  equations  are  given  and 
the  simulation  code  is  briefly  described  prior  to  a  description  of  the 
results.  Extensive  results  are  presented  which  describe  the  scaling  of 
the  fluctuation  levels  at  the  initial  stabilization  and  in  the  asymptotic, 
quasi-steady  state,  regime  with  linear  growth  rate.  In  Sec.  Ill,  we 
describe  the  pertinent  theoretical  considerations  including  a  brief  dis¬ 
cussion  of  the  linearized  growth  rate  of  the  nonlinear  interaction,  and  a 
presentation  of  the  phenomenological  model  which  describes  the  scaling 
laws.  A  summary  and  discussion  appears  in  Sec.  IV. 
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II.  The  Numerical  Simulation 


The  basic  equations  to  be  solved  are  well  known 


where  E(k)  and  6n(k)  are  the  mode  amplitudes  of  the  electron  and  ion 


oscillations,  y  is  the  growth  rate  of  mode  k  of  the  beam-driven  electron 


plasma  oscillations,  v  ^  ar>d  are  the  total  damping  rates  (collisional 
and  collisionless)  due  to  electrons  and  ions,  uu  2s  4ne2n0/m,  m  and  M  are 


the  electron  and  ion  masses,  and  c  is  the  ion  acoustic  speed.  In  the 


derivation  of  (1)  and  (2),  we  followed  Zakharov's  procedure  and  (1) 
averaged  over  the  fast  time  scale  u)g  ,  (2)  neglected  the  electron  non¬ 
linearity,  and  (5)  described  the  motion  of  the  electrons  and  ions  on  the 


basis  of  the  hydrodynamic  equations  with  phenomenological  damping 


These  equations  emphasize  the  importance  of  the 


coupling  between  the  electron  and  ion  modes.  Note  that  Eqs.  (1)  and  (2) 


are  the  k-space  representations  of  the  equations  used  in  the  study  of 


plasma  solitons 


In  the  usual  quasi linear  theory,  the  second  term  on  the  right-hand 


side  of  Eq.  (l)  is  absent,  since  only  electron  modes  are  considered.  As 


> 

a 

* 

t* 
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a  result,  stabilization  is  achieved  only  when  y^  =  0  and  a  plateau  in 
velocity  space  is  formed.  However,  consideration  of  the  effect  of  the 
ponderomotive  force  due  to  the  electron  plasma  oscillations  (i.e.,  the 
right-hand-side  of  Eq.  (2))on  the  ions  permits  the  self-consistent  ex¬ 
citation  of  ion  acoustic  oscillations.  Thus,  stabilization  can  be 
achieved,  for  y^  >  0,  when  the  level  of  ion  oscillations  reaches  a  level 
at  which  the  right-hand-side  of  (1)  vanishes  or  becomes  negative.  It 
should  be  noted  that  in  models  in  which  the  ions  are  treated  as  a 
motionless,  neutralizing  background,  the  right-hand-side  of  (2)  vanishes 
due  to  the  M-3  dependence  and  this  effect  disappears. 

Eqs.  (l)  and  (2)  describe  the  time  evolution  of  the  spatial  Fourier 

spectrum  of  the  electric  field  and  ion  density  fluctuations.  Under  the 

assumption  of  periodic  boundary  conditions,  the  spectrum  goes  from  a 

continuous  to  a  discrete  one  and  the  nonlinear  coupling  terms  become 

sums  rather  than  integrals.  Further,  because  of  the  increase  in  damping 

with  decreasing  wavelength,  the  short  wavelength  terms  can  be  neglected 

and  a  truncation  of  the  infinite  sum  over  modes  is  possible.  The  system 

of  ordinary  differential  equations  which  results  has  been  integrated 

numerically  using  a  deferred  limit  integrator  routine  with  an  adjustable 

convergence  parameter.  Numerical  tests  were  performed  to  insure  the 

insensitivity  of  the  physical  results  to  variations  in  the  numerical 

parameters.  In  general,  a  mode  spectrum  -.34  <  k\  <  .34  was  found 

to  be  adequate,  with  a  separation  of  A(kXg)  =  .02  between  the  modes. 

The  electron  damping  decrement  was  assumed  to  be  of  the  form 

v  ,  ~  k  exp (-1/2  k2\  2)  for  I kX  |>  .16  and  zero  otherwise;  and  a 

non-zero  growth  rate  (y  )  was  assumed  only  for  a  single  mode, 

which,  unless  otherwise  stated  will  correspond  to  k  1  = 

o  e 
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.02, 


and  will  be  denoted  by  y^.  The  ion  damping  decrement  was  taken  to  vary 
linearly  in  k:  and  we  use  v.,  ~  m/M  1  10.  . 

Results  of  the  simulation  which  show  the  typical  time  evolution 
of  the  beam  plasma  system  are  shown  in  Fig.  1,  in  which  we  plot  (l) 
the  total  spectral  energy  density  of  the  plasma  modes  (denoted  by  W  fc), 

(2)  the  spectral  energy  density  of  the  beam  mode  (denoted  by  Wq),  and  (3) 
the  root-mean-square  fluctuation  of  the  ion  density  for  the  case  in 
which  y  / o>e  =3  x  10-4.  As  seen  in  the  figure,  the  first  stage  of  the 
interaction  is  linear  growth  during  which  the  amplitude  of  (and,  hence, 

W  ^ )  increases  rapidly  while  the  ion  density  fluctuations  remain  at  a 
relatively  low  level.  The  turbulent  phase  sets  in  abruptly  when 
exceeds  the  nonlinear  threshold.  During  this  sudden  onset  phase,  the 
level  of  ion  density  fluctuations  increases  rapidly,  and  wave  energy  is 
transferred  from  the  beam-resonant  mode  at  a  rate  faster  than  the  linear 
growth  rate.  Stabilization  is  achieved  when  the  nonlinear  energy  transfer  rate 
is  faster  than  the  linear  growth  rate.  On  a  longer  time  scale,  an  asymptotic  state 
is  established  with  slowly  varying  levels  of  spectral  energy  density  and  ion 
density  fluctuations.  The  turbulent  nature  of  the  quasi-steady  state, 
however,  is  indicated  by  the  rapid  fluctuations  in  shown  schematically 
in  the  figure  which  are  indicative  of  a  strong  coupling  between  the  modes 
of  the  system. 

In  order  to  determine  the  scaling  of  the  spectral  energy  density  in 

the  Langmuir  oscillations  and  the  ion  density  fluctuations  with  the  linear 

growth  rate  in  both  the  sudden  onset  and  asymptotic  phases,  we  have  run 

a  series  of  simulations  witli  ■>  /j  =  ,  1  -  .  1  .  The  results  are 

o  e 

displayed  in  Figs. 
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Fig.  2  shows  the  amplitude  of  the  first  peak  in  Wo(=Wto(.)ds  a  function  of 
Y  / co  during  the  sudden  onset  phase,  for  three  choices  of  k^X^  (=0.02,  0.06,  and 
0.10),  in  which  it  is  evident  that  two  regimes  exist.  It  should  be  noted  here 
that  while  changes  in  k  X  (which  correspondsto  changes  in  v  /v,  ) 

0  6  G  D 

have  an  effect  upon  the  saturation  level  of  Wq,  there  is  little  change 

in  the  overall  scaling  properties.  We  find  that  for  low  growth  rates,  the 

peak  pump-energy  scales  as  (W  )  ^  Y  ,  while  for  higher  growth  rates 

r  r  r  o  max  o 

(w  )  'v  y  The  transition  occurs  for  y  /w  -  -0003.  The  corresponding 

peaks  in  the  ion  density  fluctuations  are  found  to  scale  as  (l5n/n0)rins'v  Y0 

over  the  range,  and  is  shown  in  Fig.  3  for  k^X^  =  0.02.  The  effect  of 

variations  in  k  X  will  be  treated  in  more  detail  later  in  this  work; 
o  e 

however,  in  the  meantime  we  consider  the  case  of  kQXe  =  .02  in  more  detail. 

Fig.  4  shows  the  level  of  ion  density  fluctuations  averaged  over 

time  during  the  asymptotic  regime  (i.e.,  the  long  time  scale  result  after 

stabilization  is  achieved)  versus  y  ,  and  we  find  that  a  result  of 

o 

<(6n/n  )  >  ~  y  is  indicated.  The  corresponding  result  for  W  in 

o  rms  'o  tot 

CO 

the  asymptotic  regime  is  given  in  Fig.  for  which  we  also  obtain  that 

<  wtot/n0Te>co  ~  yQ‘  Finally,  we  plot  the  averaged  energy  deposition 

rate  (i.e.,  the  rate  at  which  energy  is  extracted  from  the  beam)  versus 

y  in  Fig.  6,  and  a  roughly  linear  relationship  is  found  here  as  well. 

'o 

In  order  to  study  the  time  variation  of  the  field  and  density 
fluctuations  in  coordinate  space,  the  spatial  Fourier  transforms  are  inverted 
at  selected  time  intervals.  The  essential  characterist ics  of  the  co¬ 
ordinate  dependence  of  the  solution  are  shown  in  Figs.  7a-7d  in  which 
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we  plot  the  ion  density  fluctuations  and  the  magnitude  of  the  electric 
field  fluctuations  as  a  function  of  position  for  the  case  which  corresponds 
to  that  shown  in  Fig.  1.  The  normalization  chosen  is  to  the  maximum  values 
attained  by  the  ion  density  fluctuation  and  the  electric  field  fluctuation 
respectively.  Fig.  7a  depicts  the  state  of  the  field  and  density  fluctu¬ 
ations  late  in  the  linear  phase  of  the  interactions,  and  just  prior 
to  the  onset  of  the  strongly  nonlinear  regime.  The  situation  shown 
in  Fig.  7b  corresponds  to  a  slightly  later  time  just  prior  to 
that  shown  in  Fig.  1  when  the  total  spectral  energy  density  reaches  a 
maximum.  As  noted  earlier,  the  processes  of  soliton  formation  and  para¬ 
metric  stabilization  are  linked,  and  it  is  evident  that  this  figure 
describes  the  early  phases  of  both  processes.  Fig.  7c  corresponds  to 
the  point  at  which  the  soliton  reaches  maximum  amplitude.  It  should  be 
pointed  out  that  the  solitons  reach  peak  amplitude  very  quickly  once 
the  nonlinear  effects  become  important,  and  that  Fig.  7c  describes  the 
system  prior  to  the  attainment  of  a  steady  state.  In  this  phase  of  the 
interaction  the  total  spectral  energy  density  is  decreasing.  Finally, 

Fig.  7d  describes  the  system  shortly  after  the  steady  state  has  been 
achieved.  The  total  spectral  energy  density  and  soliton  amplitude 
vary  little  beyond  this  point. 

In  order  to  investigate  the  effect  of  varying  choices  of  kQAe  in 

more  detail,  we  plot  tfie  initial  peak  in  W  versus  k  for  several  choices 

o  o 

of  y  in  Fig.  8.  It  is  clear  that  while  the  scaling  between  (W  )  and 
°  o  max 

—  Q 

k  varies  with  y  ,  we  typically  find  that  (W  /n  T  )  =  (k  1  )  . 

o  o  o  o  e  max  o  e 
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III.  Theoretical  Considerations 


A  phenomenological  model  of  the  onset  phase  of  the  interaction  was 

produced  in  ref.  3  by  stressing  the  analogy  in  the  linear  regime  between 

the  formation  of  cavities  by  a  modulationally  unstable  wave  spectrum  in 

the  strong  turbulence  regime  and  the  well-known  theory  of  parametric 

instabilities.  Once  the  threshold  is  reached  (i.e.,  W  /n  T  >  k‘\  2)  the 

J  o  o  e  6 

wave  spectrum  becomes  unstable  and  forms  cavities  out  of  the  uniform 
Langmuir  turbulence.  This  involves  the  transfer  of  energy  to  shorter 
wavelength  modes,  which  permits  interaction  with  the  background  plasma. 

The  principal  purpose  of  this  section  is  to  construct  a  phenomenological 
model  of  the  nonlinear  beam-plasma  interaction  which  predicts  the  scaling 
laws  between  (l)  the  fluctuation  quantities  at  the  time  of  initial 
stabilization,  (2)the  average  values  of  the  fluctuations  during  the 
asymptotic  phase  with  the  linear  growth  rate.  This  requires  a  knowledge 
of  the  nonlinear  growth  rate  and,  as  a  result,  this  section  is  divided 
into  two  parts.  In  the  first,  we  consider  the  parametric  growth  rate  in 
an  effort  to  develop  some  insight  into  the  scaling  of  the  nonlinear  growth 
rate  on  the  wavelength  and  wavelength  spread  of  the  unstable  spectrum  as 
well  as  on  the  effect  of  electron  and  ion  damping.  Using  the  estimates 
and  scaling  found  in  the  first  part,  we  develop  the  phenomenological 
model  of  the  nonlinear  stages  of  the  interaction  in  the  second  part  of 
the  section. 

A.  The  Linearized  Dispersion  Equation 

Useful  insights  into  the  nature  of  the  nonlinear  stabilization 
mechanism  can  be  obtained  by  consideration  of  the  beam  generated  waves 


as  a  pump  spectrum  with  frequencies  in  the  vicinity  of  the  plasma  fre¬ 
quency  and  with  wavelengths  in  the  range  kQ  -  AkQ/2 <k<  k0+AkQ/2 .  We 
then  examine  the  linear  stability  of  the  spectrum  using  Eqs.  (l)  and  (2). 
This  is  analogous  to  determination  of  the  growth  rates  of  the  parametric 
instabilities  excited  by  the  oeam  modes,  and  stabilization  of  the  beam- 
plasma  instability  can  be  expected  when  the  nonlinear  term  on  the  right- 
hand-side  of  1  dominates  over  the  term  in  y^(i.e.,  the  linear  growth 
rate ) . 

The  linear  dispersion  equation  takes  the  form 


(JL'(  t 


+  iv . , 
lk 
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M  e 


f 


dk* 


W  (k*) 

o 

n  T 
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j  UJ  2(kX  )4  - 
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(3) 


where  WQ(k’ )  =  |E(k')|p/8rr  is  the  wave  energy  density  in  the  beam-driven 
mode a .  For  a  given  spectrum  W  (k'),  one  can  find  the  parametric  growth 
rate  as  well  as  the  regimes  of  importance  in  the  spectral  transfer  of 
energy.  In  this  section,  we  describe  the  spectrum  of  beam-driven  modes 
by  means  of  a  Lorentzian  distribution  of  the  form 


W  Ak 

w  (k1)  =  -2  - 2 

O  TT 


(k'-k  )p  +  A  k  2 
o  o 


(h) 
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It  is  clear  that  lim  W  (k‘ )  =  W  6(k  -k  ):  therefore,  such  a  model 
. ,  .  „  o  o'  o'’  ’ 

Ako-*0 

spectrum  is  useful  in  the  determination  of  the  effects  of  a  pump  spec¬ 
trum  with  both  finite  wavelength  and  finite  bandwidth  on  the  dispersion 
equation.  A  detailed  discussion  of  the  dispersion  equation  is  beyond 
the  scope  of  this  work,  and  we  shall  present  only  a  brief  discussion  of 
the  solutions  to  (3 )  and  (4)  in  the  text.  The  interested  reader  is 
referred  to  the  Appendix  for  further  detail. 

In  the  limit  in  which  k  «  k,  the  threshold  condition  for  the 

o 

oscillating  two  stream  instability  can  be  written  in  the  form 


~y~  >  3(kxJ‘ 

n  i  e 

e  e 


and  it  is  clear  that  the  effect  of  finite  k  is  to  reduce  the  threshold 

o 

required  for  instability.  It  can  also  be  shown  that  (u)  =  +  iy) 
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subject  to  the  requirement  that 


rOO4  ^  2  +  y2  » 


/  k 
(v3  r 


We  note  here  that  in  the  opposite  limit  (i.e.,  k  «  kQ)the  parametric 
decay  instability  is  obtained.  However,  the  threshold  condition  for 
the  decay  instability, 


6(k  X  )a  , 

n  T  e  M  o  e  ’ 

o  e 


is  typically  higher  than  that  for  the  OTS  instability,  and  we  confine 
the  discussion  in  this  paper  to  the  OTS  instability.  In  addition,  we 
observe  that  the  predominant  direction  of  energy  transfer  by  this  pro¬ 
cess  is  toward  k  >  k^  (i.e.,  toward  shorter  wavelengths). 

The  effect  of  small,  but  finite,  bandwidths  can  also  be  included 
analytically  when,  in  addition  to  (8),  we  have 


Ak 

y  »  3  -iT  ^e)S  ^e 


In  this  case  the  threshold  condition  is  of  the  form 


W  f  k  2  Ak  2 

-1A->5(kX)2  1-k—  +4—  , 

1  e  k 

n  p  K  -J 
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and  we  have  that  a> 

o 


is  given  approximately  by  { 14 )  and 


y ^ 


I"  ^ 


X 


(kx  Y 

e 


(11) 


It  is  clear  that,  in  contrast  to  the  effect  of  finite  k  ,  the  effect 

o’ 

of  finite  Ak  is  to  increase  the  instability  threshold, 
o 

The  scaling  of  the  OTS  growth  rate  with  the  spectral  energy  density 
of  the  beam-driven  modes  can  be  obtained  from  either  (7)  or  ( 1 1 )  in  the 
dipole  approximation  'i.e.,  kQ  =  AkQ  =  0).  In  the  regime  in  which 
5(kXe)2  <  W  /n  T  <  m/M,  the  maximum  growth  rate  is  given  by 


(y/u. 


w 

1  _ 0_ 

4  n  T 
o  e 


On  the  other  hand,  when  W  /n  T  JlkX  ,  m/M  we  find 
7  o  o  e  e  7 


y/we) 


max 


(12) 


(13) 


The  complete  dispersion  relation  has  also  been  studied  numerically 

for  cases  in  which  k  and  are  comparable^ to  determine  the  scaling  of  y 

with  k^  when  k^/k  'v  1  .  The  results  are  shown  in  Fig.  9  in  which  we 

plot  the  maximum  growth  rate  as  a  function  of  k^X^  for  various  choices 

of  W  /n  T  .  We  note  here  that  peak  growth  has  been  found  to  occur  for 
o  o  e 

k  ^  2k^  in  the  cases  considered.  As  seen  in  the  figure  ,  this  scaling  is 

12 


. 


dependent  upon  W  /n  T  ,  but  y  is  relatively  independent  of  k  for 
o  o  e  max  o 

k  *  &  -°2.  However,  for  higher  values  of  k  A  (0.02  <  k  A  <0.10)  the 
u  c  o  e  o  e 

scaling  ranges  between  -v  (k^)  -18  -(k^g)-"  ,  depending  upon  Wq. 

It  should  also  be  pointed  out  that  the  transition  region  for  k  A  =  02 

o  e 

corresponds  to  those  modes  with  group  velocities  comparable  to  the  ion 
sound  speed. 

Solution  of  the  dispersion  equation  in  the  case  in  which  v.,  ,  v  .  , 

ik’  ek’ 

and  y^  are  nonzero  yields  little  additional  insight  in  relation  to  the 
complexity  of  the  problem,  and  we  limit  ourselves  here  to  a  discussion 
of  the  necessary  threshold  condition  in  the  dipole  approximation. 

The  approximate  threshold  condition, 


5(kXe)2  [l  -  4^2-  +1  (kXeP 

L  k  U) 


is  obtained  by  setting  u>  =  0 .  It  is  clear  that,  as  in  the  case  of 
finite  AkQ,  the  effect  of  damping  is  to  enhance  the  instability  thresh- 
hold,  thereby,  causing  the  transfer  of  energy  by  this  mechanism  to 
be  less  favorable  and  raising  the  saturation  level  of  the  beam-driven 
modes.  We  note,  also,  that  in  the  dipole  limit,  the  resulting  dispersion 
equation  is  similar,  but  not  identical,  to  the  one  considered  by  Nishikawa. 
The  difference  arises  from  the  fact  that  the  coupling  coefficient  between 


the  ion  and  Langmuir  modes  (X  in  Nishikawa 's  notation  ),  and  which  is  the 
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coefficient  of  the  nonlinear  term  in  (2),  was  assumed  by  Nishikawa  to 
be  constant.  Here,  the  coupling  coefficient  is  k-dependent  and  depends 
on  the  frequency  mismatch  between  the  pump  spectrum  and  the  excited 
Langmuir  waves  (i.e.,  ^  k2X‘"e  ) . 

B.  The  Phenomenological  Description 

The  phenomenological  model  is  based  upon  the  linearized  theory 
of  the  parametric  interactions.  We  consider  the  wave  spectrum  to  be 
composed  of  two  groups  of  high  frequency  Langmuir  waves,  one  of  which 
is  resonant  with  the  beam  and  is  denoted  by  and  one  composed  of  the 
nonresonant  waves  denoted  by  W  ,  and  the  ion  waves  designated  by  Wg  . 

There  are  two  regimes  of  interest  in  this  work:  the  initial  stabilization 
of  the  linear  instability  (i.e.,  the  sudden  onset  of  the  turbulence 
process),  and  the  subsequently  established  asymptotic,  quasi-steady 
state.  The  energy  balance  between  these  groups  of  spectrum  is  provided 
by  both  linear  and  nonlinear  mechanisms.  In  the  case  of  the  beam-resonant 
waves,  Wq,  the  principal  source  of  energy  is  the  linear  beam-plasma  insta¬ 
bility,  and  energy  is  lost  through  the  parametric  coupling  with  other 
modes.  The  nonresonant  spectra  (i.e.,  and  W  )  can  grow  only  by  means 
of  the  nonlinear  transfer  mechanism,  and  are  subject  to  damping  by  the 
background  plasma.  We  note  here  that  linear  damping  of  the  beam- resonant 
spectrum  is  excluded  due  to  the  high  phase  velocities  of  the  waves  involved. 

The  time  evolution  of  the  system  in  the  initial  stabilization-phase 
can  be  described  by  the  following  set  of  rate  equations 


—  W  =  2V  W 
dt  o  o  o 


2y(Wo)  wx  , 


(15) 
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On  the  basis  of  this  set  of  equations,  we  can  sketch  the  approximate 

behavior  of  the  beam  plasma  system  as  a  function  of  r.  Evidently,  Wq 

increases  linearly  with  r  (which  corresponds  to  an  exponential  dependence 

on  t)  until  the  contribution  of  the  second  term  on  the  right  hand  side 

of  (17)  becomes  significant.  The  maximum  value  of  the  spectral  energy 

density  of  the  beam-driven  modes  occurs  when  t  =  t  ,  where 
J  max' 


U)  r 
e  max 


j£n 


y 

k  -2 

e 


(22) 


Substitution  of  (22)  into  (20)  yields 


(rf)-(rf)  *  4 

\  0  0  e  /initial  e 


(23) 


Since  the  logarithmic  variation  in  (23)  provides  only  a  weak  variation 

with  y  and  'W  )  we  have  that  (W  )  ~  y  in  this  parametric 

'o  1  initial’  o  max  'o  r 

regime . 

In  an  analogous  manner,  it  may  be  shown  that  in  the  regime  in  which 

m/M,  kk  )'  <  W  /n  T  the  peak  value  of  the  spectral  energy  of  the  beam- 
e  o  e  e  r  ’ 

driven  modes  is  given  approximately  by 
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As  a  result,  the  scaling  is  given  by  Wq  «  y in  this  regime.  To 
summarize  these  results,  we  have  that 

(25) 

within  the  context  of  the  dipole  approximation  in  which  electron  and 
ion  damping  are  not  included.  Of  course,  it  should  be  noted  that  since 
the  inclusion  of  phenomenological  damping  results  in  a  decrease  in  the 
growth  rate  of  the  oscillating  two  stream  instability,  the  maximum 
value  of  the  spectral  energy  density  of  the  beam  driven  modes  (as  well 
as  the  level  reached  at  saturation)  will  be  higher  than  that  predicted 
within  the  context  of  the  dipole  approximation. 

Comparison  of  (25)  with  Fig.  4  shows  that  good  agreement  with  the 
simulation  can  be  achieved  by  such  a  phenomenological  model,  and  we 
observe  that  not  only  are  these  two  regimes  in  the  scaling  of  (W  )  ax 
and  y  found,  but  the  transition  point  is  also  in  correspondence  with 
the  results  of  the  simulation. 

The  scaling  between  (W  /n  T  )  and  k  X  can  now  be  readily  ex- 
a  o  o  e  max  o  e 

plained.  The  initial  saturation  occurs  when  the  linear  (y  )  and  nonlinear 

(y^j )  energy  transfer  rates  are  in  balance  (i.e.,y  ^  y^ ) .  Therefore, 

if  we  consider  the  scaling  of  the  mode  with  the  fastest  nonlinear  growth 

for  which  (see  Fig.  9)  y,„  =  (k  X  ) ’ k  (W  /n  T  ) 2,  then  it  is  clear  that 
6  'NL  o  e  o  o  e 

stabilization  will  occur  for 


which  is  in  qualitative  agreement  with  the  results  of  the  simulation. 


y0  .*<».>“  5 


,-^>3'^,  (-1  *  (kX  ) 

’  ID  M  \m)  e 
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The  asymptotic  state  is,  perhaps,  the  least  explored  regime  in  the 

study  of  this  nonlinear  stabilization  mechanism;  however,  guided  by 

the  simulation  results,  we  will  discuss  a  phenomenological  approach  in 

this  regime  as  well.  In  the  strongly  turbulent  regime  we  must  also 

include  the  low  frequency  fluctuations  in  the  analysis.  We  define 

W  =  ^(6n,/n  )s  ,  and  attempt  to  fit  the  simulation  results  to  the 
k 

following  set  of  equations,  for 


and 


W 

a  o 


at  n  T 
o  e 


=  2V 


o  n  T 
o  e 


n  T 
o  e 


.  w,  *  W 

a  i  _  *  o 

at"  nT  V  n  T 

o  e  o  e 


V  rp  * 

l  n  T  ' 

o  e 


W  . 

— —  =  2  W 
n  T  3 

o  e 


(26) 


(27) 


(28) 


These  equations  reflect  the  difference  between  the  quasi-steady 
asymptotic  regime  and  the  linear  phase  of  the  interaction;  specifically, 
the  presence  of  finite  amplitude  ion  acoustic  waves  and  cavities  allows 
a  mechanism  for  the  dissipation  of  long  wavelength  Langmuir  waves  via 

scattering  off  density  fluctuations  (i.e.,  the  Dawson-Oberman  high  fre- 

12 .  * 

quency  resistivity  ).  This  process  is  included  in  the  term  v  W^, 


where 


and  0  is  a  factor  which  depends  upon  the  detailed  character  of  the 
spectra.  Under  the  assumption  that  the  density  and  electric  field  per¬ 
turbations  satisfy  the  pressure  balance  condition,  we  expect  “  6(kXe)2 

v  ,  -f- 

from  which  it  follows  that  v  “  60  ws  Equation  (28)  follows  from 

pressure  balance  considerations.  The  damping  of  the  W  spectrum  is 
modelled  by  the  term  in  . 

A  basic  aim  in  this  paper  is  to  determine  how  well  the  phenomenological 
description  of  Eqs.  (28)-(5o)  approximates  the  simulation  results  for  the 
stationary  state.  Under  the  assumption  of  a  steady  state,  Eq.  (28) 
implies  that 

yo~3BWsie-  (30) 

As  seen  in  Fig.  4,  such  an  approximate  relationship  holds  for  the 
range  of  the  simulations  with  B  ““  .025.  Equation  (27)  yields  v 
in  the  steady  state  regime;  therefore,  the  energy  deposition  in  the 
system  (2yQw0)  will  scale  as  2yQv0  “  27  viyQ  ti  Comparison  with  Fig.  5 

__3 

bears  out  this  scaling  for  v,  “10  oo  .  This  value  of  v,  corresponds 

1  e  1 

to  Landau  damping  for  the  mode  corresponding  to  k^Xe  =»  .21  (or  an  equiva¬ 
lent  phase  velocity  of  about  5.1  v  ).  An  important  conclusion  to  be 
drawn  from  this  scaling  is  the  linear  relationship  between  the  energy 
deposition  rate  and  y  .  Finally,  Eq.  (28)  was  compared  with  the  simu¬ 
lation  results  and  agreement  was  found  to  be  better  than  40$. 
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IV.  Summary  and  Discussion 

In  the  present  work,  we  have  subjected  a  theory  of  the  stabili¬ 
zation  of  the  resonant  beam  plasma  instability  by  a  nonlinear  coupling 
of  Langmuir  and  ion  acoustic  fluctuations  to  a  detailed  analytic  and 
numerical  study.  The  interaction  proceeds  via  the  oscillating  two 
stream-instability  which  can  be  important  if  the  spectral  energy  density 
of  the  beam-generated  spectrum  of  Langmuir  oscillations  exceeds  threshold. 
At  this  point,  both  shorter  wavelength  (with  correspondingly  lower  phase 
velocities)  Langmuir  oscillations  and  ion  acoustic  fluctuations  are 
driven  unstable,  and  the  transfer  of  energy  from  the  spectrum  of  waves 
in  resonance  with  the  beam  can  become  important.  As  a  consequence,  the 
peak  value  of  the  beam-driven  waves  (and,  hence,  stabilization  of  the 
interaction)  occurs  when  the  nonlinear  transfer  of  energy  to  the  non¬ 
resonant  Langmuir  modes  and  ion  acoustic  oscillations  balances  the 
linear  growth  of  the  resonant  waves. 

It  should  be  noted  that  effects  due  to  spontaneous  emission,  stimu¬ 
lated  scattering,  electromagnetic  mode  coupling,  particle  trapping  and 
pitch  angle  diffusion  have  been  neglected  in  the  analysis.  In  addition, 
parametric  decay  interactions ,  which  tend  to  transfer  energy  to  longer 
wavelength  modes  with  higher  phase  velocities,  were  neglected  in  the 
treatment.  The  importance  of  these  effects,  however,  must  be  determined 
for  specific  applications  of  the  theory. 

The  most  important  aspects  of  this  work  are  the  determination  of 
scaling  laws  between  the  growth  rate  of  the  linear  beam-plasma  instability 
and  the  fluctuation  levels  in  the  electric  field  and  plasma  density 
attained  at  the  initial  stabilization  and  in  the  asymptotic  state. 
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as  well  as  the  formulation  of  a  phenomenological  description  of  the 
process.  Estimates  of  the  scaling  of  the  spectral  energy  density 
of  the  beam-resonant  modes  and  the  ion  density  fluctuations  are  ob¬ 
tained  by  consideration  of  the  nonlinear  transfer  rates.  To  this  end, 
we  include  a  brief  discussion  of  the  effects  of  finite  wavelength 
and  finite  wavelength  spread  in  the  pump  spectrum  on  the  growth  rate 
of  the  OTS  instability.  A  discussion  of  the  effects  of  electron  and 
ion  damping  on  the  growth  rate  is  also  included.  On  the  basis  of  this 
study,  estimates  of  the  scaling  of  the  spectral  energy  density  of  the 
beam-resonant  spectrum  at  the  point  of  initial  stabilization  are  ob¬ 
tained.  Specifically,  it  is  found  that  the  peak  spectral  energy 
density  of  the  beam-driven  modes  is  directly  proportional  to  the 
linear  growth  rate  for  Y0/we  £  m/M,  and  is  proportional  to  the 
square  of  the  linear  growth  rate  otherwise.  This  is  found  to  be  in 
substantial  agreement  with  the  results  of  the  numerical  simulation  . 
Agreement  is  also  demonstrated  for  the  predictions  of  the  phenomenological 
model  and  the  results  of  the  simulation  in  the  asymptotic  regime.  Quanti¬ 
tative  estimates  of  the  level  at  which  peak  spectral  energy  densities 
occur  are  found,  however,  to  be  sensitive  to  the  levels  of  electron  and 
ion  damping  included  as  well  as  to  the  central  wavelength  and  wavelength 
spread  of  the  pump  spectrum.  It  was  also  demonstrated  by  the  numerical 
simulation  of  Eqs.  (1)  and  (2)  that  the  parametric  saturation  mechanism 
is  intimately  associated  with  the  process  of  soliton  formation. 

An  immediate  application  of  the  parametric  stabilization  mechanism 

is  to  electron  streams  in  the  solar  wind;  specifically  to  the  case  of 

type  III  solar  bursts.  Here,  the  electron  streams  have  been  observed 
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to  propagate  large  distances  with  small  energy  loss,  which  is  in  con¬ 
trast  to  the  predictions  of  the  quasilinear  theory.  It  has  been  shown, 
further,  that  for  parameters  typical  to  such  streams  the  nonlinear 
stabilization  mechanism  described  in  this  work  is  appropriate.  Ex¬ 
perimental  verification  of  the  action  of  this  mechanism  would  be  provided 
by  the  observation  of  solitons  in  conjunction  with  the  electron  streams. 

However,  direct  measurement  of  soliton  fields  in  the  solar  wind  is,  at 
the  present  time,  difficult  and  we  must  rely  on  indirect  methods.  In 
particular,  electromagnetic  radiation  with  u  t  2u>e  has  been  observed 
to  be  associated  with  type  III  solar  bursts.  The  intensity  of  this 
radiation  is  seen  to  scale  linearly  with  the  stream  for  weak  electron 
fluxes,  and  to  scale  with  the  square  of  the  stream  density  when  the 
electron  flux  rises  above  a  certain  level.  Since  Langmuir  solitons 
may  be  expected  to  emit  2  oje  radiation  and  since  the  linear  beam-plasma 
growth  rate  is  directly  proportional  to  the  electron  density,  an  es¬ 
timate  of  the  scaling  of  the  radiation  intensity  with  the  linear  growth 
rate  can,  in  principle,  be  obtained  which  is  based  on  the  stabilization 
mechanism  described  here.  It  has  been  shown  that  the  volume  emissivity 
of  2  u>e  radiation  from  densely  packed  Langmuir  solitons  depends  linearly 
on  Wq.  Consequently,  the  radiation  intensity  is  expected  to  scale 
linly  with  the  electron  density  for  weak  electron  fluxes  and  quadra- 
tically  for  strong  fluxes.  As  mentioned  previously,  this  is  in  accord 
with  observations. 

Another  important  consequence  of  the  present  work  is  its  application  to 

plasma  heating  by  means  of  relativistic  electron  beams.  For  long  beam  pulses, 

the  energy  deposition  during  the  initial  phase  of  the  interaction  is  small 

(i.e. ,  W  /n  T  <1)  and  it  is  suffiecient  to  consider  the  asymptotic  regime, 
o  o  e 

From  Eqs.  (26)  -  (28),  it  can  be  shown  that  the  stopping  length  for  a  relativist 
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beam  with  energy  =  6mc  per  electron  (where  6  is  the  relativistic  factor)  is 
of  the  order  of 


L 


(6-  l)mc2 
T 

e 


(1  +  6). 


(31) 


The  important  observation  which  arises  from  the  above  relation  is  the  dependence 
of  the  stopping  length  on  (v^T^)  ^ •  As  a  consequence,  since  describes  ab¬ 
sorption  by  the  thermal  plasma,  it  is  clear  that  short  coupling  lengths  will  result 
from  having  beams  interact  with  preheated  plasmas.  For  example,  if  the  end  plugs  of 
a  tandem  mirror  were  preheated  to  10  keV,  the  stopping  distance  expected  will  be 
of  the  order  of  5m  for  typical  operating  parameters  (and  ~  .001aje).  It  should 
also  be  noted  from  the  preceding  equations  that  background  plasma  heating  will  be 
governed  by 


3T 
_ e 

3t 


(32) 


which  implies  that  the  plasma  temperature  will  increase  exponentially.  When 
v  /u»  ~  .001,  the  exponentiation  time  is  of  the  order  of  20yq\  which  suggests 

that  long  pulses  will  play  a  more  effective  role  in  heating  plasmas.  These 
aspects,  as  well  as  detailed  pulse  shapes,  will  be  discussed  further  in  a 
future  publication. 
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Appendix:  The  Linearized  OTS  Dispersion  Equation 

Substitution  of  the  Lorentzian  spectrum  (A)  into  the  dispersion 
equation  (3)  yields  an  integrand  with  four  poles  in  the  complex  k' 
plane  which  are  located  at 


:'  =  k  ±  iAk,  ~  k  ± 
o  o’  2 


2  (w+  iv ek-iYk) 

1  +  - - — - — 

3  to  (kl  )? 

e  e 


PI 


Whenk(Im  -  \)  >  0>  there  are  three  poles  in  the  upper  half¬ 

plane  and  one  pole  in  the  lower  half-plane.  Thus,  if  we  integrate 
along  the  real  k'  axis  and  close  the  controur  in  the  lower  half-plane, 
then  we  pick  up  a  single  residue  which  corresponds  to  the  pole  at 
k'  =  k  -  i  AkQ.  In  the  opposite  case,  we  close  the  contour  in  the 
upper  half -plane  nd  the  only  contribution  to  the  integral  is  due  to 

the  residue  of  the  pole  at  k'  =  k  +iAk.  The  result  is 

o  o 

w(u>+iv..)-k?c  2  = 

lk  s 


and  o k  =  sgn(Im  w  +  -  Y^) .  It  is  apparent  that  the  effect  of  the 

finite  bandwidth  of  the  pump  spectrum  (i.e.,  AkQ)  appears  as  an  effec¬ 
tive  damping  mechanism;  however,  the  appearance  of  the  factor 
requires  some  physical  interpretation.  In  the  absence  of  the  pump 
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spectrum  the  coupled  ion  and  Langmuir  waves  are  damped  by  the  back¬ 
ground  plasma,  and  these  waves  can  grow  only  at  the  expense  of  energy 

from  the  beam-driven  modes.  When  o,  >0  the  effect  of  the  finite  band- 

k 

width  is  dissipative,  and  there  is  effective  enhancement  in  the  damping 
decrement  of  the  coupled  Langmuir  waves.  On  the  other  hand,  when  o^<  0 
there  is  an  effective  "negative  dissipation"  with  respect  to  the  pump 
waves.  In  both  cases,  the  effect  of  nonzero  ^ko  is  to  inhibit  the 
transfer  of  energy  to  the  coupled  acoustic  and  Langmuir  modes.  As  a 
consequence,  the  broadening  of  the  pump  spectrum  is  expected  to  result 
in  an  enhancement  in  the  levels  of  the  beam-driven  waves  when  stabili¬ 
zation  is  achieved. 

In  general,  Eqs.  (A-l)  and  (A-2)  represent  a  quarter  dispersion 
equation  with  complex  coefficients 


where 


a)4  +  a.w3  +  a  w2  +  a  w  +  a  =0, 
3  2  1  o 


a3  =  -6<kXe)(Ve>“e  +  1{2r  +  vik>  ’ 


a2  =  -  [«k2  +  k2cs2  +  2v.kr  +  6i(kXe)(koXe)me(r  +  V]  , 


a  =  6(kX  )  (k  A  )w  (k2c  2  +  V..T)  -  i(v.,«2  +  2k2c2T) 
1  e  o  e  e  s  lk  XX  k  s 


a  =  k2c2B2  -  U2w4  +  6(kA  )(k  A  )m  r  k2c2  , 
o  sk  oe  e  oee  s 
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Fig.  2  —  Graph  of  the  scaling  of  the  initial  peak  (i.e.,  at  stabilization) 
WQ  with  70  for  three  choices  of  k0Ae 


30 


Fig.  7  —  Graph  showing  the  coordinate  space  reprt  entations  of  the  electric  field  and 
plasma  density  fluctuations  (a)  in  the  linear  phase  of  the  interaction  just  prior  to  the 
onset  of  the  turbulent  phase,  (b)  when  the  peak  in  W0  is  reached, (c)  at  the  time  of  max 
imum  amplitude  of  the  spiky  turbulence,  and  (d)  during  the  asymptotic  phase  of  the 
interaction. 


Fig.  7  —  Graph  showing  the  coordinate  space  representations  of  the  electric  field  and 
plasma  density  fluctuations  (a)  in  the  linear  phase  of  the  interaction  just  prior  to  the 
onset  of  the  turbulent  phase,  (b)  when  the  peak  in  W0  is  reached,  (c)  at  the  time  of  max 
imum  amplitude  of  the  spiky  turbulence,  and  (d)  during  the  asymptotic  phase  of  the 
interaction. 


